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ABSTRACT 

We re-investigate the old problem of the survival of the five globular clusters orbiting 
the Fornax dwarf galaxy in both standard and modified Newtonian dynamics. For 
the first time in the history of the topic, we use accurate mass models for the For- 
nax dwarf, obtained through Jeans modelling of the recently published line of sight 
velocity dispersion data, and we are also not resigned to circular orbits for the glob- 
ular clusters. Previously conceived problems stem from fixing the starting distances 
of the globulars to be less than half the tidal radius. We relax this constraint since 
there is absolutely no evidence for it and show that the dark matter paradigm, with 
either cusped or cored dark matter profiles, has no trouble sustaining the orbits of 
the two least massive globular clusters for a Hubble time almost regardless of their 
initial distance from Fornax. The three most massive globulars can remain in orbit as 
long as their starting distances are marginally outside the tidal radius. The outlook 
for modified Newtonian dynamics is also not nearly as bleak as previously reported. 
Although dynamical friction inside the tidal radius is far stronger in MOND, outside 
dynamical friction is negligible due to the absence of stars. This allows highly radial 
orbits to survive, but more importantly circular orbits at distances more than 85% of 
Fornax's tidal radius to survive indefinitely. The probability of the globular clusters 
being on circular orbits at this distance compared with their current projected dis- 
tances is discussed and shown to be plausible. Finally, if we ignore the presence of the 
most massive globular (giving it a large line of sight distance) we demonstrate that 
the remaining four globulars can survive within the tidal radius for the Hubble time 
with perfectly sensible orbits. 



1 INTRODUCTION 

' The Fornax dwarf galaxy is the largest of the dwarf 
spheroidal galaxies of the Milky Way. Its integrate d stel- 
lar l uminosity is L„ = (1.58 ± 0.16) x lO'^L© (iMated 
1 19981 ). but it has retaine d virtually no mass (< SOOOMq) 
in gas l|Mateo et al.|[l99l] ). which makes pinpointing its ra- 
dial velocity less a ccurate. Regardless, the best estimates 
l|Mateo et al.lll99lD suggest a Heliocentric distance and ad- 
vancing radial velocity of 138 ± 8 kpc and 53 ± 3 km s~^ re- 
spectively. Fornax has a well measured surface density pro- 
file enabling the deduction of its King model parameters 
which define the number density of stars from the centre 
to the tidal radius, Rt — 2.08 ± 0.18 kpc. This is comple- 
mented by the concentration p arameter log^ n(Rt/Rc) = 0.72 
l|lrwin fc Hatzidimitri"oulll995^ . 

Stellar population synthesis models tell us that for the 
ages and metallicities of stars in For nax, the ratio of mass to 
luminosity should be of order unity dKrouiillooJ), however, 
the first puzzle of Fornax, as with the other dwarfs, is that 
relating the random motions of the stars to their implied 
mass from their integrated luminosity leaves an incongruity: 



the gravity GMi,{r)r~'^ is simply insufficient. There is an 
"acceleration deficit" . 

Since we assume the measurements are correct, either 
gravity is not Newtonian, allowing the acceleration to be 
boosted according to some algorithm, or there is some dark 
matter (DM) which provides a similar service. What this 
means in detail is discussed in ^and ^ nevertheless, the 
observations of For nax have recently become so detailed 
IWalker et ahl bOOT*) that more than two thousand member 
stars can be sorted into different projected radii bins giving 
the line of sight (los) velocity dispersion (VD) as a function 
of projected radius. This basically fixes the allowed mass 
density, luminous or dark, at all radii where stars exist. 

The second puzzle is far more subtle. In orbit, there are 
at least 5 well resolved globular clusters (GCs), with im- 
portant parameters given in Table 1. The GC masses range 
between 0.37 x 10^ Mq and 3.63 x IO'^Mq, the projected 
distances are from 0.24 to 1.6 kpc and their los relative ve- 
locities with respect to Fornax are at most 8.7 ± 3.6 kms~^. 

As observers, we have no knowledge of the third spatial 
dimension (depth), nor the motion in the plane of the sky. 
What we do know, is that if we assume these GCs are a 
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GC 


Mass 10^ Mq 


Dp [kpc] 


Vios [kms '-] 


Smallest Pericentre [kpc] 


Velocity at Pericentre [kms ^] 


1 


0.37 


1.60 








2 


1.82 


1.05 


-1.2 ± 4.6 


1.20 


34 


3 


3.63 


0.43 


7.1 ± 3.9 


1.45 


31 


4 


1.32 


0.24 


5.9 ± 3.4 


1.05 


36 


5 


1.78 


1.43 


8.7 ± 3.6 


1.20 


34 



Table 1. Parameters for the 5 GCs. Masses and projected distance (Dp) come from Mackey & Gilmore (2003) and los relative velocities 
(Vios) come from Mateo et al. (1991). The pericentres and velocities at pericentre correspond to Fig llOl and are the smallest pericentres 
that can survive a Hubble time in the MOND model. 



representative sample, then there is very little chance that 
all 5 have los distances from Fornax that are considerably 
greater than the maximum projected distance, although one 
might. Therefore, it is a safe bet that all 5 currently spend 
the majority of their orbits at projected distances less than 
1.6 kpc and their orbits must reflect this. 

There is more freedom in the orbital velocity, since only 
one of the three coordinates are measured (los, not the two 
dimensions in the plane of the sky) and with far poorer 
precision than the projected radius. Three out of four of the 
GCs have los velocities of between 7 and 9 kms^^ w.r.t 
Fornax, but there is no real statistical significance attached 

to this. 

Several studies llRead et all l2006l: iGoerdt et al.l l2006l: 



Sanchez- Salcedo et al.l l2006l : Strigari et al. 12006: Inoue' 
20091 ) have suggested that the mere presence of GCs near 



Fornax; rep resents a fundamen tal problem for concordance 
cosmology l|Spergel et al."2 007^ or the modified Newtonian 
dynamics of|Milgrom (1983). They claim that the timescale 
over which the orbital angular momentum of these GCs is 
drained by the background of stars and DM, through which 
they orbit, is a fraction of the age of the Universe. Despite 
this, there is no bright nucleus of stars at the centre as there 
would be if other GCs had decayed (Tremainc 1976). This 
process which bleeds the orbital angular momentum is called 
dyn amical friction (DF) and, in the classical sense, is given 
by (jBinnev fc Tremaindlioosi ) 
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where In A is the Coulomb logarithm and is taken to be 
approximately 3. Vc(r) is the circular velocity of Fornax, 
G = 4.4x10"® kms^^fcpc^^Mg^ is Newton's constant, Mac 
is the mass of the GC and a is the measured VD of For- 
nax which is roughly 11 kms~^. We have knowledge of all 
relevant parameters in the Newtonian case once a density 
profile, p(r), of luminous plus DM has been fitted to the 
losVD. 



2 DARK MATTER MODELS 



We u se the procedure of lAngus et al.l (|2008l ) and lAngud 
l|2008h to match the losVD profile of Fornax by solv- 
ing the Jeans equation assuming a cusped NFW pro- 
file (H avarro et alJ Il997l) and the s a me cored DM pro- 



file employed by l|Read et all I2OO6I : ICoerdt et all I2OO6I : 



ISanchez-Salcedo et al.ll2006l : llnouj2009l : hereafter R06, G06, 
S06 and 109 respectively). We neglect the density corre- 
sponding to the stars, since it has no direct bearing on the 
dynamics and can be simply subtracted if one wants to infer 
the exact density of DM, which is not our goal here. The 
fits to the losVD are shown in Fig [T] (solid line) and the 

NFW density profile (^p{r) = ^.^^"^"^^ ^ is plotted in Fig [2] 

with parameters — 2.5 kpc, po ~ 1.2 x 10^ MQkpc~'^ , a 
concentration of c = ^^2M — 7,0 and the velocity distribu- 

tion of stars is isotropic (/3(r) = 1 — = 0.0). The enclosed 
mass of stars and DM is given in Fig|3l 

It is this constraint, and not the mere existence of GCs 
in orbit that must fix the DM halo since the NFW pro- 
file works perfectly well in this case and it is not possible 
to have a better match to the data, regardless of whether 
7.0 is a sensible concentration parameter for this size of 
halo, nor if the DM density is suspiciously low for a cold 
DM particle. Nevertheless, to compare with the work of 
R06, G06 and 109 we take the cored density profile they 
use (p(r) = [1 + (77)^] "^ where po = lO^^A'/ofcpc"^ 
and Ts — 0.91 kpc) and attempt to match the losVD. It is 
not possible to obtain a good match even with a variable 
anisotropy and to have a satisfactory fit one needs highly 
radial orbits, /3(r) — 0.46. A more aesthetic cored profile 
has parameters po — 2.2 x W^MQkpc~^, rg = 0.5 kpc and 
isotropic velocity anisotropy (/3(r) = 0), but we only men- 
tion it here. In any case, as one sees from Fig [5] which plots 
the deceleration from dynamical friction given by Eq[TJ the 
two DM profiles are indistinguishable beyond 0.3 kpc and so 
we discuss them jointly hereafter. 

Knowledge of the mass profile also tells us the circular 
velocity profile (Vc(r)^ = J^^ p{f)-4:Tvf^df), thus, we have 
all the pertinent information that allows us to follow the 
orbits of the GCs. 

This case where the Fornax dwarf is dominated by the 
existence of DM particles, as opposed to modified gravity, 
requires very little discussion because it has absolutely no 
problem with sustaining the orbits of GCs for a Hubble time. 

Our procedure is to numerically solve the equations of 
motion for a GC with a starting x and y position and velocity 
w.r.t Fornax (we use zero z distance - line of sight - and 
velocity for ease) . The acceleration due to gravity (Fig |3| 
is computed from the circular speed (Fig[6ll, Vc{r)^r~^ , and 
the deceleration from dynamical friction is given by Eq[T]and 
shown in Fig [5l Combining all this allows us to follow the 
orbit of the GC using simple Eulerian time-stepping such 
that for the x coordinate 
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3 DYNAMICAL FRICTION IN MOND 



Vein) X, ^ , .K/ 



(2) 



where the time step dt — 0.01 Myr, is less than 1 per cent 
of the dynamical time at 1 kpc. To update the y coordinate 
we simply swap the x subscripts for y. 

Although in the DM scenario we use only circular orbits, 
the use of both x and y coordinates (and not merely r) allows 
us to use more exotic orbits like highly elliptical ones. 

In Fig[7]we show the decay of the GCs from a distance of 
1.4 kpc from Fornax which is close to the maximum observed 
projected distance. This starting distance is nothing more 
than a sensible, assumed current 3-d distance for the 5 GCs 
to illustrate a point. Clearly, GC 3 decays to the centre in 
less than 5 Gyrs, but it takes about 10 Gyrs for GCs 2 and 5, 
which is a large fraction of the Hubble time. GCs 1 and 4 can 
survive in orbit for a Hubble time if they begin at 1.4 kpc. So 
really, if the GCs began their orbits at 1.4 kpc then only one 
of the 5 GCs has a problem surviving, although actually they 
need to do slightly more than survive because the current 
projected distances are around 1 kpc from Fornax. 

Fig|S]shows orbits that start at projected distances that 
decay towards Ikpc from Fornax (the average current pro- 
jected distance). GC 3, the most massive one, must origi- 
nally fall from outside Fornax's tidal radius to currently be 
near 1 kpc, which is absolutely no problem. In any case, as 
we shall see in the MOND section, non-circular motions can 
slow the infall rate by allowing the GC to spend significant 
portions of the orbit in weak DF regions. Furthermore, since 
the other 4 GCs are not at all vulnerable to DF, the GC that 
is becomes a statistic of one and very little significance can 
be attached to it. For instance, this single GC could easily 
have a los distance of several kpc and may not even be bound 
to Fornax. Therefore, the existence of these GCs cannot be 
used to gainsay the type of DM halo surrounding Fornax: 
whether it be cored or cusped, the orbits are perfectly stable 
to DF. 

The reason it has been mistaken as a problem for the 
concordance cosmology is due to the initial conditions cho- 
sen by the three main previous studies (R06, GOG, S06) 
and m ore recently by 109. Without justification. iRead et al.l 
1I2OO6II begin all sim ulations from a starting distance of 1 kpc, 
iGoerdt et all ||2006D increase this slightly to 1.1 kpc for their 
maximum starting di stance (but often much lower) and 
ISanchez-Salcedo et al.l (2006) decrease it again to 0.9 kpc, 
but often use 0.6 kpc to emphasise their point. 109 begin 
most simulations at 0.6 kpc, and some at 1 kpc. 

Sticking to the facts, there is no evidence that the GCs 
were at a distance of Ikpc a Hubble time ago. The only 
clues are that there are now 5 GCs with a maximum pro- 
jected distance of 1.6 kpc. It is only the statistical signifi- 
cance of observing 5 within a projected distance of 1.6 kpc 
that makes us surmise that the los distances are not sig- 
nificantly larger than the projected distances. Nevertheless, 
there is nothing to say they could not have decayed from a 
larger projected distance to their current one. What these 
authors have shown conclusively is that starting with circu- 
lar orbits from less than 1 kpc is not a favourable method 
to create the current positions. 



When iMilgrornI (|l995l ') first studied the dynamics of the 
dwarf galaxies in MOND there were only single value (cen- 
tral) losVDs and therefore only a rough guide to the internal 
acceleration of each dwarf was required to predict the M/L 
because other unknowns such as the velocity anisotropy, lu- 
minosity, distance amongst the still unknown impact of in- 
terloper stars made it a rather inaccurate business. In that 
case, the dwarf was deemed to be either isolated, wherein its 
internal gravity was stronger than the external field of the 
Milky Way, or dominated by the external field of the Milky 
Way. Now that the velocity dispersion of Fornax along the 
los is measured fairly accurately as a function of radius, it is 
important to include the MW in the MOND equation since 
it is a non-linear gravity theory. 

Therefore, we must numerically solve for the gravity at 
all radii from Fornax using the equation 



g{r)ti{x) = GAL{r)r~ 
where j-i{x) — ^ 

(170 kms^^)^ 



(3) 



and 



Rmw 



a detailed, MOND-inspired Milky Way model and lXue et al.l 
2008 for the evidence suggesting the flat rotation speed of 
the Milky Way is 170 kms~^). Notice that in the absence of 
the external fleld, g^x from the Milky Way, we would regain 
the standard MOND relation, g — ^g^ao, for g « Uo. 

Mi,{r) (plotted in Fig[31) is the enclosed mass of Fornax 
deflned by the King luminosity profile with parameters men- 
tioned in the introduction using a mass to light ratio (M/L) 
of 1.4 whi ch is used to match the losVD profile in MOND 
(Fig □ and lAngusll2008h . 

The gravity g(r) is plotted in Fig U (dashed 
line) and is always considerably smaller than ao = 
3600(kms"^)^A:pc"^ = 1.2 x 10"^°ms"^. Taking this one 
more step, at all points in Fornax, the gravity of the GC 
is increased by the factor A^. Notice the distinction here 
between Fornax and the GCs gravity: Fornax is only mildly 
affected by an external fleld (that of the MW) and so we 
use Eq [3] to infer Fornax's gravity as a function of radius. 
The gravity of the GC is dominated by the external fleld of 
Fornax, g{r), and so we can make the approximation that 
the gravity of the GC is boosted (from its Newtonian value) 
by the factor -7^. 

With this in mind, recall that DF in the Newtonian 
sense is given by Eq[T] Analytic ally, this eq uation has pre- 
viously been altered for MOND (|Ciotti fc B inncv 200i) by 
noticing the G^ term must be corrected to account for the 
MOND effect at the GCs position i.e. in Eq[T]we must re- 
place G^j-^f (see ICiotti fc BinnevI I2OO4I and also 

ISanchez-Salcedo et al.ll2006l ). 

It should be recognised that this analytical description 
of dynamical friction in MOND is quite poorly motivated 
since this gravitational boost factor (-^fj)^ varies with po- 
sition inside Fornax, since g{r) varies with radius; actually 
by a factor of four from its maximum to the tidal radius. 

Furthermore, out to a distance \J'^^^^ from the GC itself, 
the gravity is fully Newtonian: this is around 20 pc for GC 
3. The DF formula used does not explicitly take these fac- 
tors into account, calculating the boost factor at the GCs 



= 210(kms"^)^fcpc~^ (see iMcGaughl I2OO8I for 
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position, instead of accounting for the boost factor along the 
orbit of the star that is scattered by the GC. Therefore, it is 
a slightly uncertain upper limit and ideally we would prefer 
N-b ody simulatio n s of G C orbits in Fornax. 

iNipoti et all (|2008| ) studied DF of bars in galaxies 
with high resolution N-bo dy simulations in MOND (see 
also iTiret fc Coinbe3 l2008l ). and found DF is signif i cantly 
stronger in that special case. ISanchez-Salcedo et al.l (|2006l ) 
investigated this exact issue we study and showed that circu- 
lar orbits beginning at 0.9 kpc (less than half the tidal radius 
of Fornax; and significantly less than 3 of the projected dis- 
tances of the GCs) are rather short lived and claimed this is 
a critical issue for MOND. As discussed above, this is not a 
rigorous enough treatment of the problem to label this as a 
problem for the MOND paradigm. We make a more defini- 
tive statement in Fig [9] for circular orbits beginning slightly 
within the tidal radius of Fornax (at 1.8 kpc). In particular, 
notice that the survival time for 4 GCs from this distance is 
greater than 8 Gyrs and from 1.9 kpc the decay time is infi- 
nite (discussed further in !j3.2p because the density of stars 
approaches zero. Nevertheless, circular orbits beginning at 
~ 1.6 kpc a Hubble time ago are most definitely problematic. 



3.1 Radial Orbits 

We have established that circular orbits within 1.6 kpc for 
the three more massive GCs have short survival time w.r.t 
the Hubble time, but are there other orbits that enter the 
tidal radius that can survive for significantly longer? 

Assuming the GCs do indeed enter within the tidal ra- 
dius of Fornax, the only scenario that can prevent them 
decaying in a period much shorter than the Hubble time 
is to have highly radial orbits that pericentre well outside 
the core radius, Rc — 0.4 kpc, but inside the tidal radius, 
Rt = 2.08fcpc. These orbits represent the most penetrating 
orbits that can survive a Hubble time. 

Given that the orbits extend out to several kpc, it is im- 
portant to show they remain within the inner Lagrange point 
in the GC, MW, Fornax system. A comprehe nsive analysis 
of Roche lobes was developed for MOND in IZhao fc Tiaiil 
(|2006t ). where they show a GC would have to be at a distance 

Dl « Rmw (C ivij^jV ) ^ ^- Using a Milky Way distance, 
Rmw , mass Mmw = 5 X 10^°Mq like before l|McGaugbl 
20081 and rei:erences therein) and C = 1 recommended by 
Zhao fc Tiaiil l|2006l ) for the deep MOND regime, this gives 
Dl 12 kpc which is significantly greater than the largest 
apocentre reached. 

The orbital parameters corresponding to these orbits 
are given in Table 1, including the pericentric distance and 
tangential velocity (angular momentum) at pericentre, but 
only for the four most massive GCs, since the lightest one 
is not problematic. Fig [10] shows the orbital decay of the 
GCs for the orbits with the lowest possible pericentre that 
can survive a Hubble time. Since we have only three out of 
six of the phase space coordinates we cannot with absolute 
certainty rule out the radial orbits presented in Fig IIUI In 
spite of this, the fact remains that the 5 GCs currently have 
projected distances less than 1.6 kpc (and the 4 most massive 
and significant GCs are within 1.4 kpc). Therefore, the orbits 
of the GCs must not only survive a Hubble time, but also 
reflect the current positions with reasonable probability. 



To help attach some statistical significance to this prob- 
lem, we force the orbits of the GCs to be inclined so that 
the angular momentum is entirely in the plane of the sky. 
This is not particularly convoluted because if we take the 
eigh t or more dwarf satellites of the Mi lky Way as examples 
see lLvnden-Belllll983l : lD'Onghia|[2008h , their orbits are not 
random, but follow two great circles around the Milky Way 
presumably because they were part of two smaller groups of 
galaxies that entered the Milky Way at late times. 

Therefore, we are looking to find orbits that have a high 
probability of being observed within 1.4 kpc of Fornax andwe 
calculate this probability simply by counting how many time 
steps in the orbit integration were spent within a projected 
radius less than 1.4 kpc, we found that the 4 largest GCs 
spend ~45% of each orbit within 1.4 kpc of the centre of 
Fornax. This corresponds to a probability of at most 0.45'' ~ 
0.04 for these orbits averaged over the Hubble time, which 
seems unfeasible. 

However, the special thing about these orbits is that 
within the last two or three Gyrs the GCs have spent all of 
their orbit inside 1.4 kpc. As can be seen in Fig[TO]the orbits 
are increasingly circularised with time and this amplifies the 
probability of being observed within 1.4 kpc. This time pe- 
riod can be further increased if the GCs were captured at 
some fraction of the Hubble time, since the orbits need not 
be as extreme. It must be stated that these orbital parame- 
ters are fine-tuned, but it is still important to be aware such 
orbits exist, since sadly we have no facts about their history. 



3.2 Circular Orbits at the Tidal Radius 

A crucial point is that if the GCs have circular orbits beyond 
the tidal radius (in fact 85% of Rt ^ 1.9 kpc) the orbits are 
permanently stable. As one can see from Fig[5]or Eq[T] the 
DF deceleration tends to zero around 2 kpc as the density of 
stars from the King model tends to zero at the tidal radius. 
Here, there are no stars to scatter, whereas in the Newtonian 
framework the omnipresent DM particles are still abundant. 
In actual fact, we have no dynamical tracers beyond the 
stellar tidal radius, so the DM might also be truncated. 

One way to estimate the probability of such a scenario 
is to think of the orbital time spent outside the maximum 
projected distance for a circular orbit at Dp^max ~ 2 kpc, 
where the circular velocity is 20 kms"'^. If we again accept 
the orbit is inclined so that all angular momentum is in the 
plane of the sky, the relationship for the projected position 
of such an orbit is simply 

Dp = 2sin(^) [fcpc] (4) 

where the angular frequency ^QQ^jy^ is VcD~l^^^. It can eas- 
ily be seen that for t > 77.5 Myr, Dp > 1.4 kpc which is the 
maximum projected distance of the 4 largest GCs. What 
we want to know is whether observing the 4 GCs within 
1.4 kpc is a contrivance if they have a circular orbit beyond 
2 kpc. Since 77.5% of the time a GC on a circular orbit at 
2 kpc is within 1.4 kpc, the probability of observing it within 
1.4 kpc is ("Y^^)- Therefore, the probability of observing all 

four within 1.4 kpc is (■y^f')^ ~ 0.36, which one could argue 
makes this a perfectly plausible scenario. 
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Figure 1. The data points are the losVD as measured by Walker 
et al. (2007). The soUd Hne is the Newtonian fit for a cusped 
NFW profile with = 2.5 kpc and po = 1.2 X lO'^ MQkpc'^ and 

isotropic orbits (5 = 1 — -^= 0.0. The dotted line takes the pa- 
rameters used in the analyses by Read et al. (2006), Goerdt et 
al. (2006) and Inoue (2009) with /3(r) = 0.46 to facilitate bet- 
ter agreement between the observations. The dashed line is the 
MOND fit using only the luminous matter with M/L=1.4, a ve- 
locity anisotropy that varies as /3 = —0.36 (^ ^'^ • "< — ^'^'^ 

_ (170 kms"'-)^ 



Figure 3. Enclosed mass profiles for cusped DM (solid), cored 
DM (dotted) and MOND (dashed). 
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Figure 2. Mass density profiles for cusped DM (solid), cored DM 
(dotted) and MOND (dashed). 



3.3 GC 3 

One final point worth mentioning is that the timing prob- 
lem is aggravated greatly by the existence of GC 3 (the most 
massive one by a factor of two). It is entirely possible that 
GC 3 orbits at a "safe" distance from Fornax, whereas the 
other four are not on orbits beyond the tidal radius but in- 
stead are gradually spiralling in. If this is the case then from 
Fig[9]we know that the next two largest GCs (2 and 5) will 
take at least 8 Gyrs to decay. Comparing that to the first 
panel of Fig 1101 to survive a Hubble time and plunge to 
within 1.2 kpc of Fornax the radial orbit must be relatively 
extreme. Recall that Fig [TO] was made to show the orbit 
which could plunge deepest into Fornax and still survive a 
Hubble time. If instead the orbit only skims the edge of For- 
nax, the radial orbit need not be as extreme. Thus, different 
orbits exist, both circular and exotic, which can survive in 
orbit of Fornax at the correct distance for a Hubble time. 
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Figure 4. Gravity as a function of distance for the cusped DM 
profile (solid), cored DM (dotted) and MOND (dashed). Recall 
that ao in these units is 3600. 



4 DISCUSSION 

The fact that 5 GCs currently orbit the Fornax dwarf, al- 
though it has no stellar nucleus, would seem to suggest that 
dynamical friction is a weak means of draining their orbital 
energy. Significantly, the DM profile necessary to match the 
losVD profile of the Fornax dwarf can easily accommodate 
the slow decay of the GCs to orbital distances resembling 
their current ones (FiglS}. 
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Figure 5. Dynamical Friction as a function of distance from For- 
nax and GC-3 for the cusped DM profile (solid), cored DM (dot- 
ted) and MOND (dashed). Recall that in theses units is 3600. 
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Figure 10. The orbital decay of highly radial orbits that enter well inside the tidal radius in MOND. The orbits decay in a Hubble time 
and spend a large fraction of the final Gyrs orbiting within the tidal radius. This scenario is discussed in depth in i|3.1l 
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Figure 6. Circular velocity profiles for cusped DM (solid), cored 
DM (dotted) and MOND (dashed). 




Figure 7. The orbital decay of the 5 GCs for the DM profile 
beginning with circular orbits at 1.4 kpc. 
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Figure 8. The orbital decay of the 5 GCs for the DM profile 
beginning with circular orbits at radii that enable sensible current 
orbital radii. 




Figure 9. The orbital decay of the 5 GCs in MOND beginning 
with circular orbits just inside the tidal radius (1.8 kpc). 



Dynamical friction inside the tidal radius is far stronger 
in MOND (see FigO- GCs beginning on circular orbits nar- 
rowly inside the tidal radius lose their orbital energy to DF 
and sink to the centre of Forna:x within a Hubble time. In 
the case of the largest GC, the orbit cannot be sustained for 
more than 5 Gyrs. On the other hand, highly radial orbits 
that enter well inside the luminous part of the dwarf can 
survive for a Hubble time, with the final 2 or 3 Gyrs spent 
almost entirely within the luminous part before eventually 
spiralling in. Furthermore, if these GCs were captured in 
the last 4 Gyrs, there is also not enough time to drain the 
angular momentum of even the most massive one. It must 
be cautioned, however, that these orbits are highly super- 
ficial and it seems highly unlikely that such extreme orbits 
exist at the expense of more pedestrian ones in every case, 
although it may be true for one. 

Outside the tidal radius (at 1.9 kpc) circular orbits can 
survive indefinitely, since there exist no stars to scatter and 
drag the GC. The only issue is the probability that the 4 
large GCs can simultaneously have no projected radii greater 
than 1.4 kpc, which we estimate to be greater than a third, 
which for a unique system is credible. A final approach is 
to ignore the largest GC which greatly relieves the timing 
problem, since the other four GCs can withstand the DF for 
substantially longer interior to the tidal radius. 

Taking into account the discussion above and the un- 
certainties of the analytical approach to DF in MOND, it is 
probably not worth attaching too much significance to the 
existence of these 5 GCs. The biggest mystery for the DM 
paradigm is why this was seen as a critical problem in the 
first place. 
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